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Abstract: In the time-space symmetric version of dynamical triangulation, a non- 

perturbative formulation of quantum Einstein gravity, numerical simulations without mat- 
ter have shown two phases, with spacetimes that are either crumpled or elongated like 
branched-polymers, with strong evidence of a first-order transition between them. These 

C^\ ' properties have generally been considered unphysical. Using previously unpublished numer- 

ical results, we give an interpretation in terms of continuum spacetimes that have constant 

^sO ■ positive and negative curvature, respectively in the 'elongated' and 'crumpled' phase. The 

o 

CO 



magnitude of the positive curvature leads naturally to average spacetimes consisting solely 
of baby-universes in a branched-polymer structure, whereas the negative curvature accom- 
modates easily a large mother universe, albeit with a crumpling singularity. Nevertheless, 
there is evidence for scaling in the crumpled phase, which we compare with the well-known 
scaling in the elongated phase. Using constraint effective-action models we analyze ex- 
isting numerical susceptibility-data of the phase transition and determine the behavior of 
the average Regge-curvature. We propose a renormalization of the Regge curvature and 
compare it to the curvature of the above continuum spacetimes, and also to the curvature 
implied by the Gauss-Bonnet theorem in the continuum. The latter involves a more benign 
multiplicative renormalization and suggests that simulations at larger volumes are needed 
to settle to order of the phase transition. 
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1. Introduction 

Causal Dynamical Triangulation (CDT) [|l|, [2|] appears to be a viable direction in the quest 
of a non-perturbative formulation of quantum Einstein gravity (QEG). Numerical simu- 
lations have shown the existence of a phase region in parameter space bounded by first 
and second order critical lines with non-trivial physical properties: an average geometry 
with semiclassical as well as fractal features and distant-dependent spectral dimensions 
H H? Hi ||]- A basic aspect of the formulation is a foliation of spacetime by space- like 
hypersurfaces - hence the name 'causal' - which facilitates a clear analytic continuation 
between real and imaginary time; for a review see [0]. Prior to CDT a formulation of 
dynamical triangulation that is symmetric between space and (imaginary) time - Symmet- 
ric Dynamical Triangulation 1 (SDT) |8|, || - was investigated vigorously. This direction 



1 Also known as Euclidean Dynamical Triangulation (EDT). 



to QEG was largely abandoned because the average background geometries unveiled by 



numerical simulation were deemed unphysical; see e.g. [10, 11 1 for reviews. A better un- 
derstanding of SDT results from a continuum viewpoint is desirable, and as a step in this 
direction we propose here a continuum interpretation of its average geometries. 

Without matter the SDT model has two phases, depending on the value of the bare 
Newton constant, a crumpled phase and an elongated phase. The spacetimes in these 
phases do not appear to be physical. In the elongated phase they have tree-like ('branched 
polymer') characteristics [jl2[ |^, [14], |l5|, 16 1 whereas those in the crumpled phase contain 
'singular structures', vertices and links connected to macroscopically large volumes in a 
single lattice step JIT], 18, 19, 20|. Strong evidence was found that the transition between 



the phases is of first order [21, 22 1. 

Arguments were put forward that a sufficient amount of matter might ameliorate 

|). With gauge fields the SDT model acquires a third 



the results |23|, |24| (see also 

phase, called crinkled |2(| 27 or smooth [28, 29]. This phase was also found without gauge 



fields upon adding a so-called 'measure term' |3C] to the action, with a new parameter 
that represents the number of gauge fields |26|, 27. 31]. Study of this new phase led to 
a confusing state of affairs: whereas some authors did not convincingly find physically 
attractive behavior [27, 31], others found evidence for a higher-than-first order transition 
1 28] with a susceptibility exponent suggestive of emerging graviton degrees of freedom 



| 29| , |25| . Recent studies of the third phase found that its spectral dimension also depends 
on the distance scale J32], p3jl . 

The continuum interpretation developed here is at scales fairly close to the UV reg- 
ulator. Such a situation is known in phenomenological particle theory, in the continuum 
but with a cutoff on the momenta, for example quark models of nuclear matter [p4|, 35] 



or meson models. At fixed cutoff, such continuum models can have first-order phase tran- 
sitions as a function of the bare parameters. With approximate SU(n) x SU(n) chiral 
symmetry a limit can be taken in which the masses of pseudo Nambu-Goldstone bosons 
(e.g. 7r, K, rj for n = 3) become arbitrarily small. However, other physical quantities such 
as the pion and kaon decay constants then still remain of order of the cutoff. In CDT and 
SDT, which a have remnant of diffeomorphism invariance consisting of permutations of 
lattice coordinates, a similar phenomenon may occur: massless gravitons emerging in the 
large lattice limit 2 , with a Plank length of order of the lattice spacing. We do not consider 
a first-order transition to be an obstruction to a continuum interpretation of SDT. 

Our interpretation is based on fitting continuum observables to lattice ones, similar to 
what is sometimes done in lattice field theory 3 . In earlier work together with B.V. de Bakker 



a continuum curvature was derived from a volume-distance correlator [41]. This curvature 
was found negative in the crumpled phase, positive in the elongated phase and changed 
sign near the transition. Examples of scaling were also given. The present work elaborates 



2 In this respect an interesting field-theory example is Z(n) lattice gauge theory. The model is completely 
discrete but for n > 6 it has a 'Coulomb phase' with massless photons [KM, B7], B8l B9l |40| . 

3 For example, the values of a scalar field propagator Gi(x, y) on a hypercubic lattice (x and y are integer 
4-tuples) may be arrayed as a decreasing function of r = \x — y\ and fitted by a continuum propagator 
G c (r) = zK\(mr)jr in a certain fitting domain; G c (r) is then a continuum approximation to Ge{x,y). 



on this, using numerical results obtained shortly after [41] that were left unpublished and 

incorporating also other work on scaling Jl2| and the structure of SDT spacetimes [ fJ2[ p3| . 

In section |2| we summarize some basic formulas of SDT. Section y recalls the Regge 

curvature Rr that appears in the action, the definition of the lattice geodesic-distance and 



the volume-distance correlator, the previous curvatures derived from these [41|, and shows 



some numerical results for larger volumes. In section 4.1 we mention some properties of 
smooth spacetimes of constant curvature which form the basis of a more elaborate observ- 
able (essentially a local proper-time metric) to be fitted to the numerical data. The result of 



such fits, shown in section 4.2, is a continuum curvature R c which changes sign at the phase 
transition. Section || discusses the relevance of R c to the baby- and mother-universes found 
in simulations of SDT. Evidence of scaling in the crumpled phase is presented in section pi 
and compared with that in the elongated phase. 

To compare R c with (Rr) we study in section |7] the constrained effective- action of 
the volume- averaged -Rr. This gives us a tool to obtain the average Regge-curvature - 
up to a constant - in the region of the phase transition from the data in ]j22[j. This 
is done in section || where also the constant is determined by comparison with data in 



1 44]. We propose additive and multiplicative renormalizations of the Regge curvature to 
enable comparison with continuum curvatures, such as R c and the curvature implied by 
the Gauss-Bonnet theorem. The multiplicative renormalization-constant following from 
the Gauss-Bonnet curvature comes out much closer to 1 than that following from R c . The 
constraint effective-action contains a renormalized gravitational constant 1/G which passes 
through zero at the transition. Conclusions are in section |9|. 

The continuum interpretation is guided by examples from regular lattices. While 
conceptually relevant, this material is delegated to appendix [A| since it is rather 'lattice 
technical'. In practical terms it leads to a simple conversion factor linking lattice and 
continuum geodesic distances. Appendix [B] contains alternative ways of fitting the volume- 
distance correlator with results that differ quantitatively (but not qualitatively) from those 
in the main text and in appendix ^| we give some details on the quartic-potential model 
used in the effective action. 

2. Symmetric dynamical triangulation 

In the (imaginary-)time-space symmetric version of dynamical triangulation, spacetime 
consists of flat equilateral four-simplices 'glued' together to form a simplicial manifold. On 
a manifold without boundary, the Einstein-Hilbert action without matter, 



s =' d ^I^Gf) (21) 



becomes translated into the form j|, |9[ 



S = -k 2 N 2 + k±N 4 , (2.2) 

_ 2tTV 2 __ Aq^4 + 10 6^2 , . 

87rGo 87tGo 



Here Go and Ao are the bare Newton and cosmological constants, N 2 is the number of 
triangles, N± the number of 4-simplices, 9 = arccos(l/4) and V{ is the volume of an i- 
simplex, 

tVTTT £ = VT0i; (2.4) 



i\V2 i 

£ is the edge length of the simplices and £ is that of the dual lattice. The formal grand- 
canonical partition function in the continuum is represented as a sum over all 4 triangu- 
lations T obeying manifold conditions, at fixed topology, usually taken to be that of the 
four-sphere S 4 [§, ||] , 



z 



V(diff) 



(2.5) 



-K4N4 



Z(k 2 ,N 4 ), 



(2.6) 

(2.7) 



Z(k 2 ,k 4 ) = J2e K * N >~ K * N * =J2 e ~ 
T N 4 

Z(k 2 ,N,) = J2 e K2N2 - 
T(Na) 

For large N4, the canonical partition function ( |2.7| ) behaves exponentially [45, 13 1, e.g. in 
the elongated phase 

Z{k 2 ,N 4 ) ~ (iY 4 )7(«2)-3 eK S( K2 )7V 4) Na ^ ^ (2.8) 

Hence, the grand-canonical partition function fl2.6| ) converges at large N4 for K4 > k\(k 2 ). 
For algorithmic reasons, simulations of the canonical partition function ( |2.7| ) were done 
with a term ex (JV4 — -/V4) 2 added to the action to allow for fluctuations around the desired 
number of simplices ./V4. It is then possible to select configurations at ./V4 = JV4 for the 
computation of averages. The phase transition occurs at a pseudo-critical point k 2 {Nj{)i 
the crumpled phase is in n 2 < k 2 and the elongated phase in k 2 > k 2 . 

The susceptibility exponent 7 that characterizes the leading correction to the exponen- 
tial behavior in ( |2.S| ) is found to have the branched-polymer value 1/2 deep in the elongated 
phase (k 2 3> k 2 ) [|l^, [D|, 0, p!q| , and to decrease to values near zero as k 2 \ n 2 [ 12 1 . In the 
crumpled phase the nature of the subleading behavior less clear [^, 46]. With additional 
matter 7 can be negative in the crinkled/smooth phase [26, 27, |3l], p8j, and its values at 
the crinkled/smooth-crumpled transition have been found to correspond to a coefficient in 



the conformal anomaly of continuum QEG [29, 25 [. 

The average of observables O is defined in the usual way 

1 



(O) 



Z[K 2 ,Ki 

for the grand canonical average, and 

(O) 



r 



K2N2-K4N4 



0(T), 



1 



Z( K2 ,N 4 ) 



T(N 4 ) 



k 2 N 2 



0(T), 



(2.9) 



(2.10) 



4 Every simplicial configuration is meant to be counted only once. Usually a weight factor 1/C(7~) - 
with C(T) the order of the automorphism group of a configuration - is inserted to avoid over-counting in 
discrete coordinates. 



for the canonical average. 

3. Volume, distance and curvature 

In this section we recall some observables for volume, distance and curvature that are 
invariant under discrete transformations of lattice coordinates. The basic definition of 
four-volume of a set of n^ of four-simplices is V = 714714, with V4 the volume (2.4) of a 
single simplex. In the canonical averages, the total volume is fixed to N4V4. 

A natural notion of geodesic distance is the length of a geodesic path through the 
interior of the simplices in the piecewise flat manifold. However, the computation of a such 
geodesic distance is complicated. In practice a convenient definition of geodesic distance 
between two four-simplices x and y has been used: dg(x, y) = the minimal number of steps, 
going from x to y from one simplex to the next, times L It is the minimal length of a 
path on the links of the dual lattice. A problem with d# is that lattice artifacts do not 
diminish at distances 3> £. This is illustrated in appendix [A] by simple lattice models in 
flat spacetime. 

The Regge curvature is concentrated on triangles. Its integral around a single triangle 
and its average over the volume are given by 



/ d 4 Xy/g~R = 2v 2 [2tt - 714(A) 
Ja 



f d x ^/g R 



47Tt> 2 f N 2 



me 

~2n 



Rr, 



(3.1) 



(3.2) 



where 714(A) is the number of simplices containing the triangle A. In ( |3.2| ) it is used that 
J^ A 714(A) = IO./V4, every four-simplex contains 10 triangles. At the phase transition (-Rr,) 
does not vanish, but it 'jumps' with a derivative {dRn)/dK2 that has a sharp peak. The 
form ( |3.2D suggests that R mixes with the cosmological constant (which' contribution to 
volume- averaged action is just a constant like the 9 term) under renormalization and that 
it needs an additive renormalization; it may also need a multiplicative renormalization. We 
return to these ideas in section M . 



In [41] an observable for curvature was proposed, based on the behavior of the volume 
V(r) of a ball of geodesic radius r from an arbitrary origin, in the continuum, 



V(r) 
V'(r) 



IT 



2vr 2 r 3 



1 



1 



Rr 2 
~36~ 

Rr 2 

Ha 



+ 0(r 4 ) 
+ 0(r 4 ) 



(3.3) 
(3.4) 



This was mimicked in SDT by assuming V(r) to be proportional to N(r), the average 
number of simplices within lattice geodesic-distance r. Its discrete derivative, using units 5 
in which the distance between the centers of two neighboring simplices is 1, 



n(r) = N(r) - N(r - 1) 



2_, "r,d e (x,y) } 1 

\ y I 



(3.5) 



These units £ — 1 will be used throughout this work, unless otherwise indicated; r in (3.4) and (3.5) 
have a different meaning and from section 4.2 onwards the first (continuum) one will be denotes by r c - 



is then 'the average volume at distance r\ The average is independent of x, which can be 
made explicit, e.g. in the probability to find the geodesic distance r, 



p(r) 



N, 



n(r) 



. 4 xy 



Jr,d e (x,y) 



^P( r ) 



1. 



(3.6) 
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Introducing an effective volume v e s, curvature 'observables' were based on the tentative 
correspondence 

v cS n(r) ++ V\r). (3.7) 

For this to work out well the quantum fluctuations should not be too large. The distribution 
( |3.6|) has long-distance tails caused by fluctuations like 'baby universes'. Examples of n{r) 
are shown in figure [j]. 

A fit n{r) ~ ar ?J + f3r 5 represent- 
ing the first two terms in the short- 
distance expansion ( |3.4|) gave a cur- 
vature Ry = —24/? /a, which was found 
positive in the elongated phase and 
negative in the crumpled phase, pass- 
ing as a function of k<i through zero 
near k\ pl[ . The left plot in figure 
@ shows the result of such fits for 
N 4 = 32000 and 64000. 

The curvature Ry depends sen- 
sitively on the r interval used in the 
fit. On the one hand it is desirable 
to use only small values of r in the fit 

to suppress the higher order terms in ( |3~4l , |3.7| ), but on the other hand one feels uneasy 
about the small r region because of lattice artifacts. In [41] this sensitivity of Ry was ex- 
ploited by restricting the fitting interval to just two subsequent values of r, writing n(r) = 
a(r)r 3 + /3(r)r 5 , n(r + l) = a(r)(r + l) 3 +/3(r)(r + l)) 5 , and R cS (r + l/2) = -24/3(r)/a(r). 
Plots of this 'running' curvature R e s(r) showed a minimum or approximate r-independence 
in the region 6 < r < 10. The right plot of figure || shows the values of R e g at these minima, 
Reff, miii) for the current larger volumes. This data is somewhat shifted compared to that 
of Ry but the general trend is similar. The pseudo-critical transition points as defined 
by the position of the peak in the node susceptibility are at k\ = 1.257(1) and 1.280(1), 
respectively for ./V4 = 32 k and 64 k p^] . 

In the next section the correspondence ([T7]) is explored further, guided by a study of 
regular lattices in flat space that allow for analytic calculations in appendix [A|. 

4. Geometry from the volume-distance relation 

4.1 Smooth spacetimes 

To analyze the average SDT spacetimes from a continuum point of view we compare the 
quantum expectation value of the volume at distance r with the form it takes in smooth 



Figure 1: n(r) for JV4 = 64000; red: crumpled phase, 
«2 = 1-26; blue: elongated phase, K2 = 1-29. 
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Figure 2: Left: curvature Ry as a function of K2 from a least-squares fit to the n(r) data in the 
region 1 < r < 11, for N4, = 32 k and 64 k. Right: values of i? c ff. m i n - Errorbars are obtained with 
the jackknife method. 



classical spacetimes. Because of the averaging in n{r) = \J2y^r,di(x,y)) we assume these 

spacetimes to be homogeneous and isotropic with constant scalar curvature R. Given the 

constant curvature, their Euler index \E is related to the volume V by the Gauss-Bonnet 

relation 

R 2 V 



XE 1927T 2 ' 

This can be checked by integrating the Gauss-Bonnet invariant, 

h, — it — ^MfjcJXT + -K K A^v-K ^ , 

XE = 7^2 / d 4 x ^fg E, 
over spacetimes in which locally 



from which follows that 



^k\)j,u .. „ \9nfj,9Xu 9kv9\i±) 



B-*. 
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(4.1) 

(4.2) 
(4.3) 

(4.4) 
(4.5) 



Since the SDT simulations used here have 5 4 topology we are especially interested in the 
case xe = 2. Locally, the metric line element can be expressed in terms of the geodesic 
distance r from an arbitrary origin, 



ds 2 

R{r) 



dr 2 + a(r) 2 dn 2 3 , 
a"{r) a'(r) 2 



6 



+ 



a(r) a(r) 2 a(r) 



(4.6) 
(4.7) 



with cftl| the line element on the unit three-sphere. Here a(r) plays the role of a (dimension- 
full) Robertson- Walker-like scale factor in which r is the imaginary propertime. With 
R(r) = R = constant (|4.7| ) becomes a differential equation for a(r), with the boundary 
conditions a(0) = and a'(0) = 1. For positive R the solution is 



r 12 

a(r) = r sin — , R{r) = —* 



?*0 



2 ' 



(4.8) 



where r$ is the curvature radius. The three-volume at distance r is 

V'(r) = a{rf I dU 3 = 2vr 2 o(r) 3 . (4.9) 

The domain of r can be extended to < r < 7rro, resulting in the S 4 geometry with xe = 2, 
the imaginary-time version of De Sitter space. Its four-volume is 



V 



TT/'O 



dr V f (r) = — 7T r 



(4.10) 



For negative R we have 



air) = vq sinh — , R 



-5- c 4 - 11 ) 




In this case the domain of r can be ex- 
tended all the way to < r < oo, result- 
ing in the hyperbolic space H which has 
infinite volume. Finite volume is also pos- 
sible, by division by a discrete transitive 
subgroup of the symmetry group 0(4, 1) 
of H A 1 47]. This leads to an infinite num- 
ber of nonisometric spaces, which are ob- 
tained by gluing the sides of 4D polytopes. 

These may have cusps, regions of finite volume extending to infinite distance. There are 
1171 nonisometric spaces with Euler index 1, 22 of them orientable, which have five cusps 
|48]. From double coverings, orientable 'gravitational instantons' can be obtained with 



Figure 3: V'(r) for flat 2D tori: symmetric 2x2 
(blue) and 2x4 (red). 



Euler index 2 [49]. The cusps can be 'Dehn-filled', resulting in many different topologies 
|50fj . Such finite- volume spacetimes are relevant to a semiclassical evaluation of the formal 
path integral in the continuum ]51| . 

Globally, the volume V'(r) can be a complicated function. Figure ^ shows examples 
for flat tori in two dimensions, R = 0, xe = 0. In finite volume (4.11) can only be valid in 
a finite interval, say < r < r+, and the maximal value r + will depend on the particular 
spacetime. For xe = 2 it has to be smaller than the value for which the volume of the 
coordinate patch, 



2vr^ 



r + 



8tt 2 



air 



r l 



2 + cosh 






sinh' 



4l+ 

2r ' 



equals the total volume (8tt /3)?q set by the Gauss-Bonnet formula (4.1): 



r + | XE= 2 < 2ro arccoshy / 3/2 ~ 1.317rQ. 



(4.12) 



(4.13) 



If r + would be equal to this maximal value, then the boundary region at r + - being no 
longer homogeneous and isotropic - would have to be singular, perhaps crumpled to a 
point. 



4.2 Local average geometry in SDT 

The simple relation (^^) between the volume at distance r and the scale factor in the 



metric of smooth homogeneous spacetimes, and the earlier made correspondence ( |3.7| ) 
between n(r) of SDT and V'(r), invites one using n{r) to derive an average metric from 
the SDT simulation results, 

a(r) « [n(r) Ueff /2 7 r 2 ] 1 / 3 . (4.14) 

A scalar curvature R(r) can then be computed by the standard formula ([4.7]). For smooth 
spacetimes the relation between a{r) and V'(r) can be valid only in a local patch (as figure 
|3] illustrates) and the region to the left of the maximum of the n(r) curves is the place to 
attempt the deduction of a metric scale factor from n(r). 



In making (4.14) precise it is helpful to turn to simple lattice models which can be 
studied analytically. We do this in appendix [A] for triangular and hypercubic lattices in 
flat spacetime. Applying the insight obtained there to SDT, we write 

r c = \(r - s), a c (r c ) = Xa(r - s), (4.15) 

r v & 1 1/3 

a{r - s) = — ^n{r) . (4.16) 

The subscript c denotes a continuum distance and scale factor, and A is a is a conversion 
factor that compensates for effects of using the lattice geodesic distance r. The shift s is of 
order of £, the lattice spacing on the dual lattice (we continue to set £ = 1). Discretization 
effects at small r are expected to be reduced somewhat by interpolating N(r) and replac- 
ing n(r) by N'(r - 1/2) (cf. flAlp , ( |A.20| ) and figure |l7|). We tried this with piece-wise 



polynomial interpolations of N(r) at r = 1, 2, . . . . Denoting these for simplicity also by 
N(r), its derivative is N'(r). Using N'(r — 1/2) turned out to make no significant differ- 
ence in the subsequent analysis, even when using high orders of interpolation (for linear 
interpolation n(r) = N'{r — 1/2)). The breaking of rotation invariance in SDT is probably 
less than, for example, a hypercubic lattice, similar to the notion that a hexagon is a better 
approximation to a circle than a square. The average over the random configurations will 
reduce the breaking further. 

For the models in appendix [A], v e g is defined by the continuum limit. Such a limit 
is not possible of the SDT results in this paper and v e g was determined by fitting the 
constant-curvature forms (4.8) and ( 4.1 1| ) to the n(r) data at r = r m j n , . . . , r max , as in 



(|A~2l| ). Then flQgj p gives a value of A, 

A = (V^cff) 1/4 , (4-17) 



where V4 is defined in ( |2.4| ). The effective volume £v e g is equal to the volume of a four- 
simplex scaled by the conversion factor A -4 . In the following we shall work mostly with 
the lattice r and a(r) and not with the continuum r c and a c (r c ), keeping in mind that 
continuum lengths involve the conversion factor A. 

We have done fits without shift, s = 0, called A-fits, and with s / 0, called B-fits. 
The idea behind the A-fit is to avoid lattice artifacts as much as possible and use only a 
minimum fitting domain r m ; n S> 1. Instead, with the B-type, the shift s is chosen such 



n(r)" 3 



that the data can be approximated reasonably well with least-squares fits at distances all 
the way down to r = 1, choosing r m ; n = 1. To see that this makes sense, consider figure 
^j, in which linearly interpolated nfr) 1 ' 3 data is extrapolated to the region where it has a 
zero point. The zeros points so are nearly independent of K2 and the plot indicates that 
using a negative shift s = so can improve the fits in the small distance region. Note that 
already from r = 3 onwards the crumpled- and elongated-phase curves deviate. Using s as 
an additional fitting parameter gives unstable results, the three parameters turn out to be 
too correlated. Therefore we fixed s = so in B-fits. 

The results of a A- and B-fits are de- 
scribed in appendix [B]. The choice of fit- 
ting domain and type of fit influences the 
resulting curvatures quantitatively rather 
strongly, but not their qualitative depen- 
dence on «2- Curvature is a sensitive ob- 
servable; the form ( [4.7D leads it to depend 
on second derivatives of n(r) through (4.16). 

In the following we describe a fitting 
method that may be called 'discrete oscu- 
lation' of type B (DOB), its fitting domain 
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Figure 4: Plot of nfr) 1 ' 3 , linearly interpolated 
_|_ and extrapolated into the region r < 1, for 
K% = 1.26 in the crumpled phase (upper) and 1.29 
(lower) in the elongated phase, N4, = 64 k. The 
zeros are at s (1.26) = -2.59, s (1.29) = -2.67. 



1. It is an improvement on R c g (section 
|3|) that uses constant-curvature forms for 
the fitting function with the shift s = sq. 
Specifically, in the elongated phase c and r$ are determined by the two equations 

c 1 / 3 sin[(r - s)/r ] = n(r) 1 / 3 , c 1 / 3 sin[(r + 1 - s/r ] = n(r + l) 1 / 3 , r = r min , (4.18) 

where c = 2i: 2 /v c q, and similarly with sin — } sinh in the crumpled phase. The resulting 
curvature ±12/rg is denoted 7 by R OS c(i" + 1/2) and these values are interpolated. It is 
gratifying that R osc (r) has stationary points: a minimum in the crumpled phase and a 
maximum in the elongated phase (the latter is absent in R e g(r)). By the 'principle of 
minimum sensitivity' the value of r is chosen to be the stationary point r stat of R osc (r): 
(d/dr)i? osc (r)| r=r . stat = 0, R osc = R OS c(r sta t), r 0oS c = ^/V2J\R^\. 

Figure || shows results of such a fit. In the left plot, the upper/lower data show 
nfr) 1 ' 3 in the crumpled/elongated phase and the fitting curves. Plotting the third root of 
n(r) enhances the short distance region. The hardly visible error bars are inherited from 
jackknife errors on n(r). With c = 2ir 2 /v e &_ determined, the metric a(r) follows from (|4.16j). 



Since a(r — s) = [n{r)/c] 1 '^ , it is just a scaled and shifted version of the data in the left 



plot, interpolated. The curvature R{r) of (4.7) and its components 6/a — 6a /a = R\ 
and —6a" /a = R2, R = R\ + R2, are shown in the right plot. The lower half applies to 
K2 = 1.26 were curves cross at r sta t — 7.83, the upper half is for K2 = 1.29 where curves 
cross at r sta t — 4.89. Evidently R(r) is not constant. Its component R\{r) has a region 
of slow variation around r s t a t, where it touches R osc /2. The double derivative component 

■> r m i n = r would result in matching value and first derivative at r. 



Interpolating n(r) and letting r max — > r n 
7 The r dependence here is not to be confused with that in (4.7) 
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n(r) 
15 - 




Figure 5: Left: results of DOB fits at ft 2 — 1.26 (upper data in blue and csinh[(r — s)/ro] in red) 
and at 1.29 (lower data and csin[(r— s)/ro]), for N4 = 64 k. Right: curvature R(r — s) (green-brown) 



and its components R\ = 6/a(r — s) 2 — 6a' (r — s) 2 / a(r - 
(blue) . The green horizontal lines represent R osc . 



s) 2 ] (red) and R 2 = — 6a"(r — s)/a(r - s) 



i?2( r ) varies much more. For constant curvature, the 6/a 2 cancels part of —6a' 2 /a 2 in R\ 
and the sum of the two is equal to R2 ■ 

Results at N4 = 64 k for DOB-fit parameters including other values of K2 (chosen from 
'non-outlying' data in figure |2|) are listed below, 



crumpled phase 


elongated 


phase 


«2 


s ro c 


K2 


s ro c 


1.255 


-2.58 7.50 0.0734 


1.282 


-2.65 11.0 0.121 


1.260 


-2.59 7.76 0.0727 


1.283 


-2.66 10.0 0.124 


1.266 


-2.60 8.16 0.0723 


1.285 


-2.66 9.78 0.125 


1.270 


-2.61 8.54 0.0727 


1.290 


-2.67 9.38 0.127 


1.277 


-2.63 9.88 0.0727 


1.300 


-2.68 9.21 0.126 



(4.19) 



The constant c = 2n 2 /v e g differs somewhat between both phases. It is the ratio of the 
surface of the unit 4D ball to the effective volume v e s that characterizes the average hyper- 
surface at fixed geodesic lattice distance r. The lattice-continuum conversion factor (cf. 
( 4.17 )) A = (-U4C/27T 2 ) 1 / 4 comes out as ~ 0.405 (0.465) in the crumpled (elongated) phase, 
and the resulting continuum curvature R c = R osc /X 2 = ±12/(ro OS cA) 2 is plotted in figure 
|6[ Remarkably, its linear interpolation goes through zero at the pseudo critical point 
«2 = 1.280(1) as determined by the peak in the node susceptibility in [B2[. 



5. Baby- and mother-universes 

In the previous section we found the average SDT spacetimes to be be locally close to 
constant-curvature spaces of negative and positive curvature, respectively in the crumpled 
and elongated phase. As noted in section |4.1| , in case of positive curvature the form of the 
metric scale factor ro c sin(r c /Vo c ) can be naturally extended to r c = 7rro c and the resulting 
spacetime is then globally similar to S A . However, its volume V c = (8-7r 2 /3)rQ C is much 
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Figure 6: Curvature R c from the DOB-fits as a function of K2 with linear interpolation. 



smaller than the total volume N4V4. Their ratio is 

V4N4 VefiN^ 



n c 



(8vrV3)r 4 c (8^/3)^ 



3iV 4 
4crQ 



(5.1) 



where we used (4-17). This ratio is about 49 for (K2, N4) = (1.29, 64000) in the elongated 
phase away from the transition. The discrepancy can be explained in the continuum 
interpretation by assuming that the total spacetime consists of small constant-curvature 
components 'glued' together. The gluing regions contribute negatively to the Euler number 
such that for the total spacetime xe = 2. If the volume of these 'necks' is small enough, 
their relative contribution to the Einstein-Hilbert action can still be small. The components 
are average-size baby universes, n c in total, and each of them consists on the average of 
A^/n c four-simplices (for the example at K2 = 1.29, A^/n c ~ 1300). 

In the negative-curvature case the scale factor ro c sinh(r c /ro c ) can be extended to the 
infinite domain < r < 00 corresponding to the hyperbolic space H 4 . The finite total 



volume N4V4 puts of course a limit on r c , by (4.12): 



>'0c 



2 + cosh — - l sinh - — — < n c 



2r 0c 



(5.2) 



where now n c just stands for the ratio ( |5.lD and does not have the interpretation of a 
number of components. Using the asymptotic form of the hyperbolic functions it follows 
that r c cannot be larger than r +c ~ (ro c /3) ln(32n c ), For K2 = 1.26 in the crumpled phase, 
n c ~ 182 and r +c ~ 2.9 ro c . However, according to the derivations in section ^4], the 
contribution to the Euler number of a hyperbolic ball with radius r +c follows from the 
Gauss-Bonnet formula $L%), which gives (— 12/rQ C ) 2 iV4i>4/1927r 2 = 2n c . To get down to 
Xe = 2 of the total space this has to be compensated in some way. Apparently this happens 
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by gluing to baby universes and by crumpling to 'singular structures'. Numerical studies 
of 'boundary volume distributions' [E3], and of branching-order distributions in 'minbu 



trees' [42], have led to a picture of SDT spacetimes in the crumpled phase that consists 
of one large component, the 'mother universe', which contains the singular vertices and 
links, connected to many small baby universes. In the elongated phase there is no mother 
component. This qualitative difference between the two phases appears to be a natural 
consequence of the sign of the average curvature: positive inducing small spheres, negative 
inducing large chunks of hyperbolic space. 

The transition between the two phases can be modeled by a 'balls in boxes' model 



1 5^, |53[ . The nature and entropy of SDT spacetimes has been investigated analytically in 



great detail in |13|, 14], [T^, [46]. 



6. Scaling 

In the elongated phase the components of ( |5,1| ) are supposed to be average-size four-spheres 
with small 'caps' taken out and glued along the caps' boundaries. Such a configuration 
can be mapped to a tree graph, in which vertices correspond to components and links 
to 'gluings'. A typical example is a branched polymer graph with average coordination 
number of the vertices not very different from 2. Indeed, in the elongated phase the SDT 



spacetimes have characteristics of a statistical ensemble of branched polymers [12, 16, 54]. 
One such characteristic is the value of the entropy exponent 7 ~ 1/2, another is the scaling 
behavior of n(r) fl2|| . 

In |^l[ we investigated scaling of the probability of geodesic distance r between two 
simplices, p(r) = n(r)/iV4. Suppose that a length r* can be chosen to depend on K2 and 
N4, and K2 to depend on N4, such that the following limiting procedure makes sense: 



N 4 



p* (r/r*, ^2,^4) —> /j*(a;,r), x r = r/r* — > x, N4 — >■ 00, (6.1) 



r°° 1 

/ dxpJx,T) = lim } Axp*(x r ,K2, N4) = ) p(r) = 1, Ax = — . (6.2) 
JO N^oo^ ' *-f n 

The label r distinguishes different shapes of p as a function of x resulting from different ways 
of limit taking. The scale r* was taken to be r m , the value of r where n(r) is maximal, and 
different sequences ^2(^4) were envisioned that produce scaling functions p m (x, ^2,-^4) of 
different shapes, e.g. N^j = 8000 x 2- ? , j = 0, 1, . . ., r = «2(-^4,o)- Another possibility for r* 
is the average value r* = r av = ^ r p(r) r. Assuming (x) m is finite in the limit ./V4 — > 00, the 
two are proportional for large N4: using r* = r m we have r a _ v /r m — >■ J dx p m (x, r) x = (x) m . 
For simplicity the labels K2, N4 of p and r of p will be dropped in the following. 

Figure [7| shows an example in the crumpled phase. The scaled data match within 
statistical errors and for clarity interpolated data is shown without error bars. Data were 
chosen from an available set to give approximate matching by eye. Also quantitatively they 
corresponded to smallest differences in the norm J dx \p(x) — p'(x)\, and these differences 
converge towards zero (for p m they turned out slightly smaller than for p av ). The DOB-fit 
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Figure 7: Left: Crumpled phase n(r) data at (JV 4 , k 2 ) = (8000,1.17) (blue), (16000, 1.21) (red), 
(32000, 1.23) (brown) and (64000, 1.26) (green). Right: same data scaled: r m n(r m x)/N4 = p m (x). 



parameters of the sequence are 



{N 4 K 2 ) 
(8000, 1.17) 
(16000, 1.21) 
(32000, 1.23) 
(64000, 1.26) 



S Tq C 

-2.23 11.7 0.127 
-2.25 9.90 0.112 
-2.25 8.26 0.0831 
-2.59 7.76 0.0727 



(6.3) 



The dependence of r av or r m on N4 can be fitted with a power law, e.g. 



ociV 



1/4 



with a scaling dimension d s ~ 6.2, but it can be fitted slightly better by a logarithm, 
r m = a + blnN^, which corresponds to d s — > 00; see figure 0. An infinite scaling dimension 



was also considered to be most likely in [55, 12], based on a comparison of data at fixed 
K2 deeper in the crumpled phase. The ratio r m /r av ~ 0.81, 0.82, 0.80, 0.80, is practically 
constant compared to the change in r av itself, which suggests that < x > m is finite indeed 
in the limit N± 



00. 




Figure 8: Left: r m versus ln(A^), logarithmic fit (red) r m = —10.74 + 2.803 ln7V4, and power fit 
(blue) r m = exp[1.228 + (1/6.195) ln^], for the scaling sequence in figure 0. Right: Curvature 
radii (blue) and fit (red) via tq = r m XQ with r m taken from the logarithmic fit in the left plot and 
x a = 1.49/(-7.16 + ln7V 4 ). 

The evidence for scaling presented in jD]] was compatible with the idea that sequences 
^2(^4) could be found for which the curvature radius ro grows with r m , and that p would 
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have a semiclassical form for all x, with d s — > 4. The subsequent discovery of the first order 
nature of the transition between the crumpled and elongated phase made such an outcome 
unlikely, and instead of growing one expects the tq to approach a constant in lattice units. 
The curvature radii of the DOB-fits in table [T^ are decreasing, but they can be fitted well 
by a function that approaches a finite limit 8 (~ 4) as N4 — > 00 (right plot in figure |8|) . 

The near constancy of tq can be used to argue that negative curvature is the reason 
for the logarithmic dependence of r av and r m on N4. Since the mother universe has a 
macroscopic fraction of the total volume its average local curvature will be close to the total 
average - negative with radius ro - and it will dominate the shape and properties of the 
distribution n{r) in regions with substantial probability. One such property is the position 
r m of the maximum of n{r). Consider the integral p = J^ , dx p m (x) = (I/./V4) J, m dr n(r), 
which is almost independent of N4 in the scaling sequence: p ~ 0.34247, 0.34248, 0.34249, 
0.34240. We write p as a fraction / of the integral obtained by replacing n(r) by its 
constant-curvature fit in the whole region r < r m : 



f f rrn r — r" 1 3 
p = — / dr c ro sinh 

^4 J\ [ r 



f Cr V e 3(r m -s)/r 



24iV 4 



1 + 0( e - 2 (^-s)Ao) . ( 6 .4) 



Then 

tq ro 24p 



ni _ lnAT 4 + _^i n - +s . ( 6 .5) 

•J d J CTq 

Assuming that / is depends only modestly on N4, less than linear, the leading dependence 
is given by the explicit lnA^. Its coefficient ro/3 is for the scaling sequence 3.9, 3.3, 2.8, 
2.6, not far 9 from the fitted coefficient 2.8 of figure ||. 

Unfortunately, similar matching data in the elongated phase is inaccessible to us now. 
In [|12j, r&v/N^' was found to become independent of K2 deep in the elongated phase and 



7-1/2 



scaling was observed with r* = Nj , with a scaling function p[x) corresponding to generic 
branched polymers (no need for a parameter r). The scaling function can be characterized 
by a so-called Hausdorff dimension dn that is identified from the small x behavior 

p{x) ocx^ 1 , x^0. (6.6) 



General arguments based on a scaling assumption [12] relate the large distance behavior 
to the scaling dimension of r* oc N 1 ' , 

p(x) oc x a exp f — cix^- 1 J , x — >■ 00, (6.7) 

and it is assumed that dn = d$. For branched polymers the natural definition of distance 
is the number of links between two vertices, and an analog PBp( r ) ~ the probability of 
distance r between two vertices - scales in the generic case with d s = dn = 2 and a = 1 



1 16, 54]. With a suitable choice of r* we may write, 

p{x) = xe~ x ' , generic branched-polymers. (6-8) 



For the B-fit with c fixed by the result in table B.2 the ro are already almost constant (cf. table 
The A-fit deteriorates for the two smaller values of -/V4. 
9 For the B-fit the numbers are 3.0, 3.0, 2.9, 2.9. 
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Evidence was presented in [jig ] that this form also applies to the n(r) of SDT sufficiently 
deep in the elongated phase and good fits of ( |6"7^ ) to n(r) were obtained for all r>l, with 

1/2 

a = 1 and d s = 2. The fits used only one parameter, c\, with r* = 7V 4 , since the nor- 
malization of p is fixed by J dx p(x) = 1. Rewriting this as n(r) = (N^/rf) r exp(— r 2 j2r\ ), 
the parameter t\ = r* , which is also the position of the maximum of this analytic form for 
n(r), and it would be fixed to T\ = r m if the fit were perfect. We find that such a 'zero 




Figure 9: Left: Volume-distance correlator n(r) (blue) in the elongated phase ( (A^,K2) = 
(64000, 1.29)) and fit ( |6.9| ) (magenta). Right: close-up showing also the asymptote (iV-i/r^) (r — si) 
(dark grey) and the S 4 DOB-fit c(r sin[(r — so)/ro]) 3 (red). 



parameter fit' is poor to the data at (N^,K2) = (64000, 1.29) (the same values as used in 
figure 7 in Jl^]), and even keeping n free does not give very good fits. However, with a 



small shift in the fitting function the linear behavior (6.6) can be made compatible with 
the data, i.e. 

r — s\ 



nyr) = — xe 

n 



x 2 /2 



n 



n. 



(6.9) 



In this way we obtained good fits in the region r > 7 to the data at «2 = 1-29 and 
1.3. The first case is shown figure ||, for which s± ~ 2.2, r\ — 68 (r m ~ 73). The slope 
at r = s\ and the DOB-fit are also shown in the right plot, which exhibits a smooth 
transition from curved to linear behavior. Closer to the transition, K2 < 1.285, the form 
p(x) = xexp(— x 2 /2) is not able to fit the data well anymore even with the shifted variable 
x = (r — s\)/r\. Presumably, the crumpled phase is too near for these K2 and the two-state 
nature of the first order transition makes itself felt by contaminating the statistics through 
crumpled-like configurations. 

In the elongated phase the components in ( |5.1[ ) are small average-size baby universes 
which corresponds a vertices of a branched-polymer graph. The relation between the SDT 
lattice-geodesic distance and the branched-polymer distance is somewhat fuzzy, but the 
same limiting scaling function appears to emerge. In the crumpled phase there is in addition 
a large mother-universe component; how this reflects on the p observable is unclear. As 



mentioned in section ^Jj, there are many possibilities of finite- volume hyperbolic spacetimes 
and an analytic scaling form may not exist or be very complicated. In the following we 
compare qualitative features of p in both phases. 
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Figure 10 shows scaling functions at k<i = 
1.29 and k 2 = 1.26 for iV 4 = 64000, using r* = 
r m to match the maxima at x = 1. In the elon- 
gated phase the 'one-component region' r < r$ 



corresponds to x < xq = ro/r„ 



0.128. 

The inflexion point in x < 1, the point of max- 
imum slope, is not much larger, x~ s = 0.144. 
Presumably xq and x~ s vanish like l/r m oc 
oo 



r l/2 



12] and there is no rea- 




Figure 10: Scaling function p m (x) for «2 = 
1.26 in the crumpled phase (blue) and 1.29 in 
the elongated phase (red), N4 = 64 k. 



l/iV 4 1/z as N 4 

son to doubt that the pure generic branched- 
polymer form p m (x) = xexp(— x 2 /2) emerges 
in the limit. 

In the crumpled phase the scaled curva- 
ture radius is larger, xq = 0.38. For the other three members of the scaling sequence in 
figure [/] the scaled radii are xq = 0.82, 0.60, 0.45 (increasing N4). Since we expect the r$ 
to approach a finite limit as N4 — > 00, the xq will vanish like 1/V m , i.e. only logarithmically 
oc 1/ In N4. They can indeed be approximated by the function xq = 6/(c+ln N4), a fit gives 
(6, c) = (1.47, —7.16). This gives also a fit to the ro of the scaling sequence using only 
two parameters via ro = r m xo, which is undistinguishable from the right plot of figure || 
The inflexion points are much larger than in the elongated phase, for the scaling sequence 11 
x~ s ~ 0.6628(10), 0.6776(3), 0.6966(2) and 0.7025(1). They seem to approach a finite limit 
and can be fitted well by the form a + b/(c + In N4). However, the three parameters are 
too correlated to get a meaningful result if we want to extrapolate beyond N4 = 64000. 
Assuming a = 1 gives a good fit with (b, c) = (—4.96, 5.77). There is also an inflexion point 
x+ s in the region x > 1; its scaling sequence 1.342(21), 1.334(25), 1.308(49), 1.333(9), will 
be discussed below. 



i*p> 



ln(p) 




Figure 11: Left: \np m (x) as a function of x for k^ = 1.26 with exponential fit in 1. 
(red) and Gaussian fit in 3.5 < x < 4.9 (blue-brown). Right: same as a function of x 2 . 



< x < 3.0 



Turning to the long-distance region x 3> 1, the form (6/7) with d s — > 00 implies 



exponential behavior [12|. The left plot in figure 11 shows lnp m , which can be fitted well 
by a function linear in x in the interval 1.8 < x < 3.0. For larger x there appears to be 
a turnover to a different behavior. The right plot shows ln/> m versus x 2 , which suggests 

10 DOB-fit ro values, unless otherwise mentioned. 

n The error bars correspond to the envelope of the Jackknife errors of the n(r) data. 
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Gaussian behavior at large x and a good fit lnp = &o — x 2 /2x 2 can be obtained in x 2 > 12, 
or x > 3.5. For comparison both fits are shown in the left and right plots. The turn-over 
from exponential to Gaussian behavior is so 'abrupt' that a linear + quadratic fit to \ap 
in the whole region 1.8 < x < 5 does not look convincing. We saw similar behavior of In/) 
for the other members of the scaling sequence. The fitted values of x\ are 1.29, 1.00, 0.94, 
1.05, with errors that are hard to quantify as they are dominated by systematics related 
to the fitting domain. 




Figure 12: Left: p m (x) for (k 2j N 4 ) = (1, 16000). Right: hxp m (blue) with fit b - x 2 /2x\ in 
2.1 < x 2 < 6.6 (blue-brown), as a function of x 2 . 

Accepting tentatively the Gaussian behavior of the tail of the distribution, its inter- 
pretation might be that even in the crumpled phase, the baby universes are sufficiently 
many and small in size to cause branched polymer-like features at the largest distance 
scale. This idea can be tested by going deeper in the crumpled phase, at smaller curvature 
radii. The best example in our data set involves somewhat unfortunately a smaller number 
of simplices, (n2, N4) = (1, 16000). Since k^ (16000) ~ 1.22, K2 = 1 is here much deeper 
in the crumpled phase than the data we have shown thus far. Its DOB-fit parameters are 
s = —1.98, c = 0.143 (A = 0.360), vq = 6.03, and the maximum of n(r) is at r m = 13.7. De- 
spite the relatively small volume the putative number of components ( |5.1j ) is fairly large 12 , 
n c ~ 63. Its scaling function p m is shown in figure [t^; xq = 0.441, x~ ns = 0.7846(6) and 
x ms = 1-210(2). The right plot shows Gaussian behavior, ln/9 m is linear in x 2 over eleven 
e-folds in the tail region x 2 > 2.1 (x > 1.5). We found no convincing 13 indication of 
exponential behavior in x > 1. 

Consider again the right plot in figure [/]. The scaling violation in x < 1 shows a 
systematic steepening of the scaling function with crossings near x = 0.8 and an increasing 
maximum, with a probability f^, dx p m (x) that stays nearly constant as N4 in creases, as 
noted above. The systematics in the region x > 1 is somewhat less clear but the curves for 
the largest two N4 cross again near x = 1.6, indicating a narrowing of the distribution. If 
this trend continues the scaling function might start looking like the one in the left plot of 
figure 12. The true limit ./V4 — > 00 might even be a Dirac delta function (note the difference 



in hight of the maxima in figures [7] and [j^) 



12 The B-fit gives even n c ~ 126. 

13 A linear function of x fits the data well only in the rather small region 2.14 < x 2 < 3.34 (1.5 < X < 1. 
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To investigate this quantitatively we ex- 
trapolated the characteristics x^ s and x\ 
to larger volumes to see whether they ap- 
proach those of the scaling function in fig- 
ure [l2|: pm . The x~ s of the scaling se- 
quence have the smallest errors. They can 
be fitted by x" s = 1 - 6.11/(9.54 + lniV 4 ), 
and extrapolating to larger JV4 this matches 
the x~ s value of p m ' at In ^4 = 18.8 
(N 4 w 1.5 x 10 8 ). The fit to the x data 
does not share such a matching property, 
but xq characterizes scaling- violation. In 
x > 1, fits to the x+ s and x\ data do ap- 
proach those of pm shifted to In ^4 = 
18.8 but they miss by quite a lot, see fig- 
ure 14 [13|. Hence, the three extrapolations in this figure do not appear to support the idea 
that a unique scaling function emerges in the crumpled phase at very large N4. 

However, the reason might be finite size effects. The first order nature of the phase 
transition was seen only in the data at N4 > 32000 and N4 = 16000 may be too small to 
give a good indication of the limiting scaling function, no matter how deep in the crumpled 
phase. The finite-size effects appear to be larger in the region x > 1 than in x < 1, since 
in the former we could match x~ s at 'a reasonable' lnA^ = 18.8, whereas matching ie+ s 
would require 'an unreasonable' \11N4 ~ 71. Furthermore, we found no evidence of an 
exponential region in the tail of p m ' , instead it is Gaussian. For p\a ' ' the tail is 

mostly exponential 15 . 

In practice, the crumpled-phase scaling function is not unique because of the logarith- 



Figure 13: Fits to the scaling sequence of £+ s 
(blue, upper), x~ s (blue, lower), x\ (brown, mid- 
dle) , and extrapolations to the region around the 
shifted data of (k 2 , N 4 ) = (1, 16000). 



mic slowness of its change with N4, and the shape-distinguishing parameter r in (6.1) is 
needed indeed. 



7. Constraint effective action 

We wish to compare R c with the average Regge curvature (-Rr,) = (Rr) (the volume- 
averaged Regge curvature Rr was defined in ( |3.2|) ) . An overview of the latter in a relatively 
large region of K2 values can be found in (5^, |30|, Q|, but here we need a closeup near «§• 
It can be obtained by integrating the susceptibility ((Rr — (Rr)) 2 } oc (d/dK2)(RR.), for 
which suitable data at N$ = 64000 is shown in [^] . We shall do this by fitting a model for 
the constraint effective potential of Rr to a set of data values. 



14 Because of the substantial systematic uncertainty in the Gaussian fits leading to the Xi data points, 
these have not been given error bars but large dots, and the fit to their scaling sequence is of the least-squares 
type. 

15 The associated probabilities are small: /~ dx p£' 16000) (x) ~ 0.024; for (« 2 , N 4 ) = (1.26, 64000) the 
probability in the exponential region 1.5 < x < 3.5 is ~ 0.127, in the Gaussian tail x > 3.5 it is ~ 0.0014. 
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The distribution of R R is described by the constraint effective action T, 

/oo 
dRe~ r W, (7.1) 

7 (1V4) "°° 

where the delta function stands for applying a Dirac function after interpolating a histogram 
of R R . Properties of the constraint effective action in scalar field theory and its relation 



with the usual effective potential are discussed in |57|]; see e.g. |>g] for a numerical study 



and [59[ for a study of a composite gauge-invariant operator (<p*<p) in electroweak theory. 



Rewriting ( |3.2|) in the form 

iV 2 = ^— (r r + 209-) , V = N A v A , (7.2) 

4TTV2 V V ±J 

we see that e K2 2 can be taken out of the SDT sum, such that we can write 

T = 5 + S, S = k r V(-R + 2A ), KR = -^ = —^—, Ao = -lO0^, (7.3) 

47TU2 167rGo V4 

E = -In J2 S(£ 2 R R ,PR). (7.4) 

T(JV 4 ) 

Here, S is the Einstein-Hilbert action specialized to constant curvature R. Note that E, 
the constraint selfenergy function, does not depend on k. It follows from these definitions 
that 

R, v = (Rr) = ^-^-lnZ + 2A , (7.5) 

V dn R 

X^V((R 2 R )-(R R ) 2 ) = ±-^lnZ. (7.6) 

One expects E/V to depend only moderately on the volume. For large V the integral in 
(7.1) can then be done in a saddle-point approximation. For a first order phase transition E 
is supposed to have two minima, not necessarily at equal depth. Although E is independent 
of k, the minima in T do depend on k; at a certain value k r they are at equal depth, i.e. 

KR (-R + 2A ) + ^^ = 1 + ^-(R-R ± ) 2 + O((R-R ± f), R^R ± , (7.7) 

where we assume + and — to correspond to the elongated and crumpled phase, respectively. 
We then have 

'" - ( R - R , Y R x< R , i: = ;,-,, - /,•;. - 

16vrG 

(7.8) 

Here k can be seen to define a renormalized large-scale G. Near k = we have to keep 
both saddle points, 

Z « e ~y^^+i) Zi z = z++z _ } (7.9) 

/ ~4\ 1/2 

I 2tts ± £ \ e v iKR±+K 2 s±/2 ) _ ^_ 10 ^ 



r(i?) w F 



«(--R + 2A ) + 7 + tt 1 - (^ " fi ±) 2 

ZS;)- 



*±=/? j dRe V[nR-{R-R±f/^) = 



V J 

p(R) = P+ {R)+p-(R), p±(R) = I e V[^-(H-H±) 2 /2 S ±] ) (7 . n 
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where p±(R) are the probabilities of R in the elongated (+) and crumpled (-) phase. Using 
the notation 

R± = Rs±R d , s± = s s ±s d , 8= -la.—, w = kR a + K 2 s d /2, (7.12) 

4 s_ 



the average curvature and susceptibility following from (|7.q ), ( |7.q ) and (|7.9| ) can be ex- 
pressed as 

i? av = R s + ks s + (i?d + KSd) tanh(u;V + 5), (7-13) 

X = X P + Xb, (7.14) 

X P = F(i? d + Ks d ) 2 [1 - t&nl 2 (wV + 5)] , (7.15) 

Xb = s s + s d tanh(wy + <J)]. (7.16) 

The susceptibility \ splits naturally into a peak component Xp that is proportional to 
the explicit V and a background Xb that does not have this V dependence; the other 
parameters are expected to become volume-independent for large V . The value k c where 
X P is maximal is close to zero, 



6 I 3 1 -2 \ 1 ~/,,_ 3 x 



+ (z't-^ 2 i> d )- s - ? ;+0(V- 3 ). (7.17) 



With x given, i? av in (7.13) is its first integral with respect to n, in which R s plays the role 



of integration constant. 

It turns out that the Gaussian model ( fflj| ) is not quite good enough to fit the sus- 
ceptibility data outside the peak region. An extension, in which the generating functions 
In z± are obtained from a Legendre transformation of an effective potential that includes 
also quartic terms, is able to give a good fit with only two new parameters t±. At kr = k r 
the potential near the minima is 

f±(R) = ^ {r ~ r±)2 + S^ {r ~ r±)4 - (r - i8) 

Solving the 'semi-classical equation' [d/dR)[nR — f±{R)] = gives R = R± ln (K), and 

lnz± = Vw±(k), w±(k) = KR s f n (K) - f±(R!f n (K)). (7.19) 



For t± = this reduces to the Gaussian model without the root prefactor in ( 7.1Q 1 ), such 



that 5 = in ( |7.13| ) - ( 7.17 ). Details are in appendix ^|. The variance w±(k) falls like 



k 2 ' 3 for Kt± S> 1, unlike the Gaussian model for which it is constant, w±(k) = s±. The 
series expansion in k, 

W± ( K ) = K R ± + k 2 s± Q - 1 kHI + A K H% + 0(k 6 4)) , (7.20) 

shows the behavior at small k. The difference of w± with that of Gaussian model is only 
a function of nt± . 
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8. Phase transition and renormalized Regge curvature 



In this section we renormalize the average Regge-curvature (Rr) (|3.2| ) and compare it with 
the continuum curvature R c of figure |6| and with the continuum curvature that follows from 
the Gauss-Bonnet formula fl4.1| ). This Gauss-Bonnet curvature, Rqb = ±(3847r 2 /V) 1 ' 2 , is 
that of a smooth spacetime of constant curvature and volume V = N4V4. 

The Regge curvature can be compared with 'the plaquette' in SU(N) lattice gauge 
theory. Classically, the Wilson loop of the discretized gauge field around an elementary 
square (x,fj,u), TrU^ ux , is related to the field strength F^ u (x) appearing in the continuum 
action, by 4^ Tr(l — U^ ux ) = £ 4 F fll/ (x)F fJil ,(x) + 0(£ 6 ). In the quantum theory, the 
numbers change: for an expectation value in the ground state at zero or finite tempera- 
ture, 4^Tr(l - (U^ x )) = Z + Z F 2£ A (F^{x)F^ v {x)) + 0(£ 6 ), where Z and Z F 2 are 
dimensionless renormalization constants which depend on the gauge coupling and the tem- 
perature. In the standard model there are in addition to the gluon condensate (F 2 ) also 
condensates bilinear in the quark fields, and the Higgs condensate. We assume a similar 
property for the Regge curvature, 

R, v = (Rr) = z r 2 + z R R rcn + o{F), (8.1) 

where R re n is a renormalized 'Regge condensate'. The 0{£ 2 ) is somewhat inappropriate, 
since a continuum limit is not available at this stage. In QCD renormalization 'constants' 
like Zq and Zpi are computed in perturbation theory, which serves to define the non- 
perturbative condensates (see for example |K| and references therein; a similar example in 
electroweak theory is in |^9[). Here we shall turn this around and consider -R ren as given, 
and see what this implies for Zr, after having defined Z$. Our candidates for R Ten are R c 
and i?GB- First we need to obtain R av from the susceptibility data at the phase transition. 



Indications of a first-order nature of the transition were presented in [21, 22]. A two- 
state signal was found in the average number of nodes (Nq), for JV4 = 32000 and 64000. 
Furthermore, the exponents A and T characterizing the growth rate and width of the peak 
of its susceptibility 

2 \ / AT- \ 2 \ 



«-*<#>-<£> 



XN oc N£, 5k 2 oc iV 4 " r , were estimated to be A = 0.81(4) and T = 1.24(18) (22|. For a 
first-order transition these exponents are 1. 

The number of nodes at fixed ./V4 is equivalent to the volume-averaged Regge curvature 
( |3.2j ) because of the relation iVo = N2/2 — N4 + xe, with xe = 2 for spherical topology. 
Hence, the node susceptibility is proportional to the variance of the volume-averaged Regge 
curvature, i.e. the susceptibility defined in (|7.6|), 



XN = c xx x, c xx = - ^ ~ 1.97 x 10 4 . (8.3) 



The left plot in figure 14 shows a least-squares fit of the quartic model ( 7,18| , T.1S| ) 



to x^o data at JV4 = 32000 and 64000, taken from figure 6 in [22|, with the constraint 
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Figure 14: Left: Fit of the Gaussian (dashed) and quartic (fully drawn) model- % to \n data 
(blue) from p2]; red: c xx x, magenta: c xx Xb, low peak: N4 = 32 k, high peak: N4, — 64 k. Right: 
corresponding subtracted Regge curvature i? av — R s for N 4 — 64 k. 



that the ^-parameters are the same for the two volumes. For comparison, also shown are 
Gaussian- model fits (dashed), for which x becomes constant away from the transition. 
This quartic-model fit to the combined data looks good enough. The constraint on the 
^-parameters reduces correlations between the s- and t-parameters in the fit. Without 
the constraint, the resulting fits look 'perfect', but the parameters s+ and t + differ wildly 
between the two volumes. The fitted parameters depend moderately on the volume: 

£ 2 R S £ 2 R d s+ t 2 t + s_ I~H- k* 2 4 

0.113 288 717 1.25905 1.25712 

0.107 375 4308 798 1574 1.25898 1.25683 

0.117 299 864 1.28038 1.27976 
7.06 0.116 366 4308 887 1574 1.28029 1.27971 



7V 4 
32000 
32000 
64000 
64000 



•4) 



The value of R s was obtained from figure 2.3 in [44], which shows that (N2/N4) ~ 2.41 at 
K2 = 1.3 and N4 = 8000, . . . , 32000. Assuming the same value for 64000 this gives i? av ~ 
7.28 £~ 2 with (3^), and since the quartic fit has given the combination [i? av — R s ]k 2 =i.3, 
we get R s ~ 7.06 £~ 2 . The deviations ±i?d from this central value are smaller by almost 
two orders of magnitude. The dimensionless distance 2vYfid = \V{R+ — R-) between 



the minima of the critical effective potential, Vf±(R) = ^ f R\JV — R±\/V 
becomes larger than the largest width, J~sZ w 30, only for ./V4 > 8000. 



l^/V-R±W) /4 + 



The background curves Xb of the two models are also shown in figure |14[ Especially 
for the smaller volume the Gaussian-model background is substantially smaller at the 
transition than the quartic one. On the other hand, the fits to the total x appear equally 
good in the transition region. Hence, systematic errors are larger in Xb than in %. The 
background-subtracted peak of the susceptibility, x ~ Xb, is equal to R\V at k = 0, i.e. 
K 2 = K 2' which is close but not equal to the pseudo-critical point k 2 , the position of 
the maximum where dx/dn2 = 0, which is listed also in ( |S.4p . After subtracting the 
background, the 'improved' exponents A come out larger than one 

r 1 peak 

[X — Xb] 64000 



ix 



-ipeak 
XbJ 320OO 



2.11, 2.29, A~ 1.08, 1.19, 



(8.5) 
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respectively for the Gaussian and quartic model. Similar ratios at k\ follow from the ratios 
of R 2 d V: A~ 1.11 and 1.24. 



The right plot in figure 14 shows the Regge curvatures corresponding to the left plot 
for N4 = 64000, shifted by Rg, i.e. i? av — R s which vanishes at k = 0. For the discussion 
of the renormalization to R c it is useful and interesting to make a fit also to the R c data. 
The Gaussian model can be used for this purpose. The form (7.13) for R av is by itself not 
able to fit the -R c data well because its slope away from the transition (where tanh — > ±1) 
is tightly coupled to the width of the transition. However, with an additional overall scale 
factor 1/Z' , i.e. 

R c = (1/Z') [R' s + s' s k' + (R' d + s' d n') tanh[F(i^«' + s' a k' 2 /2]] , k = (k 2 - 4*)/(4tto 2 ), 

(8.6) 
a good fit is obtained, as shown in the left plot of figure ^jj. For comparison we have also 
shown R av . Away from the transition the Ht-R c becomes linear in K2 and the same holds 
for Gaussian model- R av . These linear forms, 



^lin = (M Z> ) K ± K + « ± *d)«] » <,lin = Rs±R d + (S S ± 8 d ) 



(8.7) 



extrapolated to K2 = 1.28, are also shown in the plot. In the large volume limit the 
width of the transition region shrinks to zero and the Gaussian model should become exact 
(assuming the parameters stabilize in the limit). At finite volume we can then replace the 
curves towards the transition by the above linear forms. 

The renormalization constants Zq and Zr can be specified as follows. We define 
Zq = R s . Then R av vanishes at k = 0, or K2 = k\ which is almost equal to k\- Given R c , 
the factor Zr is given by the ratio (R av — R s )/R c , which is shown in the right plot of figure 
15. Towards the transition the ratio is not meaningful for a first order transition, but we can 
replace it by the ratio of the linear extrapolations ( |8.7[ ), which is also shown in the plot. The 
values extrapolated from the elongated(crumpled) side are Zr = ±Rd/((R' s ± R' d )/Z') ~ 
0.20 (0.14). 



0.5 





Figure 15: Left: Fit (blue) to the R c data (blue dots) and linear extrapolations to «;§ (black). 
Also shown is R av — R s from the right plot in figure [l| and its linear extrapolations (black). Right: 
k-2 dependence of Zr = (i? av — R s )/R c (same color scheme as in figure |lj). Also shown are the 
ratios of the linear forms (pTjj) (black). 
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The Gauss-Bonnet curvature is given by Rgb = ±12/L 2 , V = (87r 2 /3)L 4 , which 
gives I 2 Rgb ^ 0.23, 0.16 (L/l ~ 7.295, 8.675), respectively for N 4 = 32000, 64000. For 
the quartic-model the extrapolated i? av — R s = ±R<\ then give at the transition Zr = 
±Rd/RGB — 0.47 and 0.73, respectively for N4 = 32 k and 64 k; for the Gaussian-model 
these number are 0.50 and 0.73. The value 0.73 is much closer to 1 than the 0.20 (0.14) 
found above for R c case. In other words, after a multiplicative renormalization, which 
approaches 1 going from 32 k to 64 k, the subtracted Regge curvature R av — R s is that of 
a smooth spacetime of volume V = N4V4 and curvature Rqb- 

9. Conclusions 

The continuum curvature R c is obtained from the volume-distance correlation-function by 
making a small-distance, constant-curvature, approximation which smoothes out 'distor- 
tions' caused by the lattice, and subsequently making a change of scale from the lattice- 
to the continuum-distance. The resulting R c is quite sensitive to changes in the way the 
numerical data are fitted, but the qualitative behavior of R c near the phase transition is 
robust. In this sense a qualitatively consistent picture has emerged in which curvature 
radii are in continuum units not much larger than the lattice spacing 16 . 

The positive curvature in the elongated phase leads naturally to the picture of an av- 
erage spacetime consisting solely of baby universes glued together into a branched-polymer 
structure, whereas the negative curvature in the crumpled phase accommodates easily a 
mother-universe. In such hyperbolic spacetimes of large volume, conflicts with the Gauss- 
Bonnet theorem have to be avoided somehow by singularities. This gives a new look on the 
occurrence of the 'singular' vertices and links observed in this phase. Remarkably, these 
'singular structures' seem innocuous to observables on the dual lattice and the scaling of 
the volume-distance correlator. 

The precise nature of scaling in the crumpled phase is hard to establish, because the 
slow logarithmic change of scale requires exponentially large volumes. Surprising is the 
branched-polymer character of its scaling function at the largest distances 17 . 

The phase transition data of J2^] could be well described by models of the constraint 
effective-action, which also gave analytic expressions for the background of the susceptibil- 
ity peak. The model with Gaussian potentials in the two phases appeared sufficient within 
the peak region and more robust than its extension to quartic potentials, which was needed 
for a good description outside the peak. Subtracting the background from the peak led 
to critical exponents closer (and even larger) than the 1 of a first order transition 18 . The 
renormalized coupling 1/G in the effective action changes sign at the transition. Of course, 
this G need not be the same as the renormalized Newton constant Gn characterizing the 
strength of the gravitational interaction. A simple adaptation also gave a good fit to R c . 
Its rapid passing through zero as a function of 1/Gq is caused by the mixed contribution 
of both phases, away from the transition it is slowly varying. 



16 Recall^ = y/Wt~3.2t 

17 The associated probability is small. 



18 



Since we refrained from estimating errors all numbers should be viewed with caution. 
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The quartic-model fit led to a detailed description of the average Regge curvature R a 



through the transition, up to a constant R s which could be found from [44]. Subtracting 
this constant, a renormalized Regge curvature R ren = (i? av — R s )/Zr emerges that passes 
through zero at the transition. This was compared with R c , and also with the Gauss-Bonnet 
curvature -Rgb , the continuum curvature that relates volume and Euler index by the Gauss- 
Bonnet theorem. The multiplicative renormalization constant Zr in this comparison came 
out rather small compared to 1 for R c at the transition, but much closer to 1 for Rqb, and 
it increased with volume. 

The average Regge curvature is quite different from R c in its sensitivity to quantum 
fluctuations. It can be expressed as a volume average of the curvature at a triangle and in 
this sense it is associated to a large distance-scale, whereas R c is derived from an averaged 
quantity at a scale not much larger than lattice spacing. The Gauss-Bonnet curvature is 
associated with a length scale L derived from the volume V = N4V4: L = [V r /(87r 2 /3)] 1 ' 4 , 
Rgb = ±12/L 2 . It is in a sense a 'target' curvature that one would like to get out of a 
detailed understanding of the renormalization of the curvature scalar in the bare Einstein- 
Hilbert action. 

The fact that the subtracted i? av — R s is not very different from Rgb for the largest 
simulation example analyzed here (N4 = 64 k) is intriguing and results at larger volumes 
are needed to establish that it is not a mere coincidence. If Zr would keep approaching 
the target 1 , then the phase transition would not be of first order 19 . 

In any case, a better understanding is needed of the physics of the branched-polymer 
structures, which are also present in the spatial slices of the spacetimes in the 'De Sitter 
phase' of CDT §. 
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A. Simple lattice models in flat spacetimes 

In this appendix we study the volume-distance relation for a 2D triangular lattice and 
for hyper-cubic lattices in D = 2, 3 and 4 four dimensions, and also the retrieval of the 
metric scale factor a(r) of the line element ds 2 = dr 2 + a{r) 2 d^l 2 D _ 1 . In flat spacetime, 
a(r) = r and this is to be the target result in the limit of zero lattice spacing, the continuum 
limit. Ensembles of these lattice configurations are supposed to contain only one member, 
so there are no effects due to averaging. 



Figure 16 shows a triangular lattice in which the lattice geodesic-distance from an 
arbitrary origin is given by the numbers in the centers of the triangles. In these units, the 
link length t = y/?>. The centers at even r form a hexagonal shape; at odd r the shape 
is somewhat different, but this becomes negligible as r increases. However, the difference 
between the lattice geodesic-distance r and the continuum distance d c depends on direction 
and does not diminish as r — > 00. For example, along a link direction, d c = -^ r (even r), 

19 Assuming that the jump 2i?d of _R av at the transition stabilizes in lattice units, as expected for a 1st- 
order transition, a doubling of the volume would bring Zb, within a few percent of 1, a further doubling to 
about \/2, etc.; in the contrary case R^ oc 1/L 2 for large L. 
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whereas in a direction perpendicular to a link direction (a direction along a dual-lattice 
link), d c = | r for r = 4, 8, 12, ... and d c = | r + ^ for r = 1, 5, 9, .... Inspection shows 
that n(r) is given by n(r) = 3r, r = 1, 2, . . . . With N(0) = 1 it follows that 



r o o 

N(r) = 1 + ^ n(r') = 1 + -r + - r 2 , r = 0, 1, . . . 



(A.l) 



As polynomials in r, n{r) and A^(r) are naturally extended from integers to real numbers. 

To avoid confusion we indicate a con- 
tinuum object by a subscript c, for exam- 
ple r c , V c (r c ) = itr 2 . Lattice objects could 
be given the subscript £, but for notational 
convenience we drop this, so r; = r. Con- 
sider the volume (area) v 2 N(r). Making ex- 
plicit the lattice spacing £ and the spacing 
of the dual lattice I = £/y/3, the dimension- 
full lattice geodesic-distance r is an integer 
multiple of £ and the volume of a triangle 

v^3 ^2 



IS V 2 



IV3£ 2 . Then 



v 2 N(r/£) 



in the limit 



4 V 2 



+ -ri+i 2 




(A.2) 



0, r fixed. It has the 



Figure 16: Equilateral triangular lattice with 
lattice-geodesic distances indicated up to r = 5. 



same behavior, oc r 2 , as V c oc r 2 . Equating 
V c (r c ) = lim^ t>2 N(r/£), it follows that 



Xr with vrA 2 = § v 2 /£ 2 , A = (9V3/8TT) 1 / 2 ~ 0.79. As might be expected, this factor A 



lies between the direction coefficients | 



0.75 and ^ ~ 0.87 found above. 



At finite £ -C r, v 2 N(r/£) ~ Vc(Ar). This approximation can be improved by including 
a shift of order of the lattice spacing, chosen to cancel the leading 0{£) contribution. This 
shift turns out to be 1/2 in lattice units: 



£ 2 N{{r/l) - 1/2) 



r 2 + -t 



(A.3) 



In fact, not only the leading lattice artifact linear in £ is canceled, but also the coefficient 
of the £ 2 term has been reduced by the shift. Furthermore, for a slowly varying function 
F(r), F(r) - F(r -£) = l(d/dr)F(r - \I) + 0(£ 3 ), and using this for F(r) = N(r/l) we 
have 



n(r/I) = N(r/£) - N{{r/£) - 1) = 
= N'({r/I)-l/2) + 0(£ 3 ) 



l± N{(r/£) - 1/2) + 0(I 3 ) 
dr 



For the simple triangular case there are no 0(£ ?y ) corrections: £N'((r/£) — 1/2) 
ln(rjl). 



(A.4) 
3r = 
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Using r c = Xr and V c (r c ) ~ V2N({r/£) — 1/2), the metric scale-factor a c (r c ) can be 
retrieved from n or N' through 

^-V c {r c ) « A-^NHr/l) - 1/2) = ^|iV'((r/£) - 1/2) = ^n(r)(A.5) 
dr c Adr \£ \£ 

(A.6) 



27ra c (r c ) 



27rr c . 



In this 2D example there are no finite £ corrections, a c (r c ) = r c . 

The constant A can be absorbed by changing units. Suppose the function V c , which 
in D spacetime dimensions has engineering dimension D, is described in terms of a length 
scale roc, e.g. a curvature radius. We can then introduce objects without the subscript c, 
V and a, that 'look' exactly like their continuum version V c and a c , 

V c (r c ;r 0c ) = X D V c (r;r ) = X D V(r,r ), r = r c /X, r = r 0c /A, (A.7) 

a c (r c ;r 0c ) = Ao(r;r ). (A. 8) 

In short, V c (r c ) = X D V(r), a c (r c ) = Xa(r), and we have to keep in mind that distances 
derived from r are larger by a factor 1/A than those corresponding to r c . In terms of a, 
the generalization of ( |A.5| ) to a lattice in D dimensions is tentatively 

DCnair) - 1 = Dc D X^ D ~^ a^r.) - 1 = X-^-^V^r.) (A.9) 

(A.10) 

(A.H) 



VD 



£X D 

C D = vr, 4vr/3, vr 2 /2 



N'{{r/£) - 1/2), 

D = 2, 3, 4, 



where Cd is the volume of the unit ball in D dimensions. The factor A is determined 
by the coefficient a that characterizes the small-distance behavior of N after taking the 
continuum limit (cf. the discussion after ( A.2| )), 

ID ATt(„ ie\ i /n\ „D , <nt„D+2\ f A 12) 



lim £ u N{{r/£) - 1/2) 

v D N((r/£) - 1/2) 

A 



ar~ +0(r u+z ), 



V c (r c 

( VD 



D+2\ 

D„D | /nf„D+2 7>2\ 



(\ 



C D X u r u + 0(r 

l/D 



(A.13) 
, , (A. 14) 

\e D c D J 

The order of the corrections 0(r +2 ) and 0(£ 2 ) applies to the lattices in this appendix. 

The tentative equations (|A.1C| )-( [A,14| ) above can now be turned around to define a 
metric scale factor a(r) from the lattice N(r) at finite lattice spacing. For reasons to 
become clear below and in the main part of this paper, we shall allow for a shift of order 
of the lattice spacing in the relation between r c and r, 

r c = X(r -s), s = 0(£). (A.15) 

Then 



a c (r c ) = Xa(r - s), 



a{r 



VcS 



VcS 

Dc D 

VD 

£X D ' 



N'((r/i) - 1/2) 



1/(D-1) 



or 



VcS 

Dc D 



n{r) 



l/{D-l) 



(A.16) 

(A.17) 
(A.18) 
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is a concrete realization of ( 4,14 ). Note that 



Hr-s)]^ = r + O(r 3 ), (A.19) 

in accordance with the proper-time interpretation of r in the continuum. The version with 
n(r) in ( A.17| ) will have somewhat larger lattice artifacts at non-zero £. 



Let us now see how this works out for a flat hyper-cubic lattice. In two dimensions, 
the centers of the plaquettes at dg = r form a square rotated by 45°. With points on 
the dual lattice labeled by integers (xi, ■ ■ ■ ,xd), the lattice geodesic distance is dg(x,y) = 
(\x\ — j/i | + • • • + \xd — Vd\)£ (£ = £)■ As for the triangular lattice, dg(x, y) differs from the 
continuum distance through the interior of the lattice, d c {x,y) = \/Yli( x i ~ Vi) 2 ^> even 
at arbitrarily large distances and there can be many lattice geodesic paths with the same 
r = dj>(0, y). Along a dual link direction, d c = r, whereas in a direction along a link of the 
original lattice, d c = r/^/2. The breaking of rotational symmetry of the set of centers at 
distance r is larger here than for the triangular lattice. For n and iV we find, reverting for 
simplicity to lattice units £ = 1, 



(A.20) 



where we used vd = £ = £ for the cubic lattices to calculate the A. Similar to the 
triangular case, n and N are naturally extended as polynomials from integers to real 
numbers. The same result is obtained by using interpolation of sufficiently high order. 
The powers of r in n{r) differ by 2, so the lattice artifacts start at order £ 2 . In contrast, 
N(r) contains all powers up to D. But for N(r — 1/2) the powers of r differ by 2 again, 
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compatible with (A.4). In D = 3 and 4 the coefficients of the non- leading powers of r in 
N'(r — 1/2) are slightly smaller than those in n(r), so the lattice artifacts in N'(r — 1/2) 
are somewhat smaller than in n(r). 



The formulas ( A.10| )-( [A.14| ) result in a scale factor a(r) with lattice artifacts in the 



small r region. These artifacts cause a non-zero scalar curvature R(r) (by ( |4.7| )) that 
rapidly vanishes as r — > oo (the continuum limit in lattice units). In four dimensions, 
choosing the version with n(r) in (| A .17 ) with shift s = 0, we have a(r) = (2r + r 3 ) 1 ' 3 , with 



27r 2 /v e ff = 8/3. The scalar curvature has the large-r expansion 20 R(r) = 8r 6 + 0(r 8 ). It 
becomes accurate for r > 3, where R(r) < 0.01, whereas R(l) — 0.74. In three dimensions, 
a(r) = (1/2 + r 2 ) 1 / 2 , R(r) = -4a" /a - 2a' 2 /a 2 + 2/a 2 = -r~ 4 + r" 6 + 0(r~ 8 ), and 
R(l) — —0.44 but already smaller than 0.004 for r > 4. It is no surprise that R(r) is of 
order of the discretization scale for r of order 1, but it is comforting to find it to be already 
less than percent of that for r > 4. In two dimensions we have simply a(r) = r, as for the 
triangular lattice. 

In 4D simplicial gravity we do not have the luxury of determining A and v e s from the 
behavior of n(r) or N'(r — 1/2) at arbitrarily large r, i.e. the continuum limit, because 



3 In the version using N'(r - 1/2), a(r) = [(7/4)r + r 3 ] 1/3 , R(r) = (49/8)r^ 6 + C(r~ 8 ). 
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the distances at which a description in terms of a metric scale factor a(r) might apply 
is not very much larger than the lattice spacing. Typically, we find in the main body of 
this paper that the maximum r has to be less than about 12. An important quantity is 
the effective volume v e g in the relation a(r — s) = [n(r)v e fi/2'ir 2 ] 1 > 3 (s is the optional shift 
introduced in ( A.18[ )). We can estimate it from intermediate distances 21 K r < r m by 
fitting a constant-curvature model for V'(r) to n(r) or N'(r — 1/2), e.g. 

3 r. 2 



n(r) 



c tq sin ■ 



To 



or 



[N'(r - 1/2)] 



2ir 2 
v e s 

1/3 , 



mini ' mm 



+ V 



1/3 ■ r 
c ' ro sm — 



ro 



(A.21) 



(A.22) 




and its sinh analogues. The sine and hyperbolic sine cases can be combined into an explicit 
function of y = i/g, with y > 0(< 0) corresponding to positive (negative) curvature. 

It is instructive to test this on the cu- n(r) w 

bic lattice case in four dimensions. We shall 
do the fitting to n(r) 1 ' 3 as in the main text. 
Figure (0) shows N'(r - 1/2) 1 / 3 , ra(r) 1 / 3 , 
their continuum limit (8/3) 1 ' 3 r, and a func- 
tion /(r) that is the linear interpolation of 
n(r), r = 2, 3, . . . , and its linear extrapo- 
lation into the region r < 2, where it has a 
zero at sq ~ —0.49, f(so) = 0. Its slope at 
the zero point, f'(so) ~ 1.27, smaller than 
(8/3) 1 / 3 ~ 1.39. Note that the curve for 
N'(r — 1/2) 1 ' 3 lies slightly below that of their continuum limit (8/3) 1 / 3 r (brown), and the 
n(r), closer to the continuum limit. Figure linear inter/extrapolation f(r) (green), for the 4D 
|j~8| show the result of two least-squares fits, CUD1C lattice, 
one with s = to the n(r) 1 ' 3 data at r = 6, 

7, . . . , 10, and one with s = sq to the data at r = 1, 2, . . . , respectively the A-fit and B-fit. 
The A-fit gives y > 0, positive curvature with a c larger than the exact value 8/3 ~ 2.67 
of the continuum limit, the B-fit y < 0, negative curvature with c smaller than 8/3. The 
B-fit shows nice agreement to the smaller-r data all the way down to r = 1, but the values 
for ro and c are actually less accurate than for the A-fit. Of course, the curvature radii 
increase to infinity and A to the value in table A.20| (c — > 8/3), as r max — > oo. 



Figure 17: N'(r - 1/2) 1 / 3 (blue), n(r) 1 / 3 (red), 



Having determined v e g we can then construct the scale factor a{r) by (A. 17) and 
compute the metric curvature R(r — s) from (|4.7j ). The latter is shown in the right plot 



of figure 18. This curvature is not constant but becomes reasonably close to zero in the 
fitting region. The plot also shows the more accurate R(r) obtained from the exact value 
c = 8/3, for which a(r) = (2r + r 3 ) 1 ' 3 , R(r) = 8r~ 6 + 0(r~ 8 ). The metric curvature R(r) 
is sensitive to the magnitude of a(r) through the term 6/a 2 in ( |4.7| ), an error in c leads to 
an error in R(r), with new terms in its asymptotic expansion starting already at order r~ 2 . 
In fact, the — 6a" /a and —6a' 2 /a 2 + 6/a 2 contributions are separately rather large but of 
opposite sign and they nearly cancel for the exact value of c. 
1 Recall that r m is the position of the maximum in n(r). 
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Figure 18: Left: two fits (red) to the n 1 / 3 cubic lattice data (blue dots); lower curve: A-fit (s = 0), 
c = 2.82 (A = 0.61), rO = 36.7 (R = +12/r% = 0.0089), fitted data at r = 7, 8, . . . , 10; upper curve: 
B-fit (s = so), c = 2.11 (A = 0.57), r = 21.7 (R = -12/rg = -0.025), fitted data at r = 1, 2, . . . , 
10. Right: curvature R(r) from the A-fit (upper red), from the exact c = 8/3 (s = 0, middle blue 
curve), and R(r — s) from the B-fit (lower red). The straight dashed lines represent ±12/r0 2 of the 
A- and B-fit. 



In this cubic case the curvature R(r) resulting from the A- and B-fit becomes accurate 
at the percent level of t~ 2 at continuum distances r c > 3.1 and 5.6, respectively. The 
accuracy of the B-fit lessens when we determine sq by linear extrapolation from = 1, 2 
instead of r = 2, 3 above. It leads to a 42% larger |so| and a slope /'(so) deviating further 
from the continuum limit. Fortunately, in the SDT case this difference is only about 15%. 

The DOB-fit plays tricks here: R sc( r + 1/2) has an apparent minimum at the left 
boundary of the r-values for which ( 4.18|) has a solution, which is r = 3. For r = 2 there is 
no solution 22 and the region r < 3 cannot be reached by interpolation, which leads one to 
reject the apparent minimum. As might be expected, the principle of minimum sensitivity 
then leads the DOB-fit to 'slide' to a stationary point which lies at infinity, r sta t = oo, 
where the fit becomes exact. For r + 1/2 = 8, the resulting R c is already more accurate 
than the one from the B-fit with fitting domain {7, 8, 9, 10}. 



B. Results of A- and B-fits 

We start with the A-fits. The left plot in figure |T^ shows the result of an A-fit of 
c[sinh(r/ro)] 3 to the n{r) data in the crumpled phase at k^ = 1.255, 1.260, 1.266, 1.270, 
for A r 4 = 64000. To avoid cluttering only the curve for K2 = 1.26 is shown in black. The 
right plot shows a result with sinh — > sin in the elongated phase (black curve, K2 = 1-29). 
In this case the fit was done simultaneously to the data at K2 = 1-282, 1.283, 1.285, 1.290, 
1.300. The parameters of the fits are c = 2ir 2 /v e s (independent of K2 within a phase) and 



2 This happens also in SDT cases, but then there are true minima in r > 3. 



31 



n(r) 
1400 




Figure 19: Left: results of fits to the n(r) data in the crumpled phase for K2 = 1.26; black curve: 
A-fit in 7 < r < 12, red curve: B-fit in 1 < r < 12. Right: similar for the elongated phase for 
k-2 = 1.29; A-fit 7 < r < 10, B-fit 1 < r < 10. Blue: numerical data with jackknife error bars. 

the ro depending on k 2 ; their values are in the following table: 



crumpled phase 




elongated phase 


c=0.305 


(A = 


.436) 


c=0.323 


(A = ( 


K 2 


ro 




«2 


ro 


1.255 


9.54 




1.282 


8.28 


1.260 


10.4 




1.283 


7.43 


1.266 


11.9 




1.285 


7.21 


1.270 


13.5 




1.290 


6.94 


1.277 


40.7 




1.300 


6.85 



0.442) 



A-fit, 



0. 



(B.l) 



The lattice-continuum conversion factor A corresponding to c is also listed in (B.l). The 



deviations in figure 19 at r < 6 are larger in the elongated phase than in the crumpled 
phase (but note that the vertical scale is enlarged by an order of magnitude) . 

In this respect, the B-fits that include the shift s do visually much better. Figure 
2C shows results of a B-fit with fixed s = —2.63, which is the average of the two so 
of figure |3]. The same k 2 an d N4 were chosen as in figure [l^. In the crumpled phase, 
c i/3 ro s inh[(r — s)/ro] was fitted with the least-squares method 23 to ra(r) 1 ' 3 ; in the elongated 
phase sinh — > sin. For comparison we have also shown the result of the A-fit in figure 20 
for k 2 = 1.255 and 1.300 (dashed curves), and on the n(r) plot in figure 19 the result of 
the B-fit (red curves). 

The B-fit gives clearly a better description of the data. (However, for the 4D cubic test 
case in appendix [A| the A-fit is more accurate than the B-fit in estimating v e g and in repro- 
ducing the zero curvature.) It is comforting to note in passing that two-parameter ('devil's 
advocate') fits of the form ar + /3r 3 to n(r) 1 ' 3 at individual k 2 give larger least-squares 
deviations than two-parameter (c and ro) B-fits. Note also that the B-fits reasonably ap- 
proximate n{r) over a larger domain than the osculation fits in figure H The B-fit analog 



23 0nce s is fixed by the condition n(s) = 0, similar results can be obtained with chi-squared fits to n(r) 
data in the region r > 7. 
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Figure 20: B-Fit to the data in the crumpled phase (upper set of red curves, fitting domain 
1 < r < 12) and elongated phase (lower set of red curves, fitting domain 1 < r < 10), with fixed 
so = —2.63. The black dashed curves represent the A-fit for K2 — 1.255 and 1.3 . 



of the A-fit table in (glj) is 



crumpled phase 




elongated phase 


c=0.0853 


(A = 


.317) 


c=0.115 


(A = 0.341) 


K 2 


''0 




K 2 


TQ 


1.255 


8.31 




1.282 


13.9 


1.260 


8.65 




1.283 


11.6 


1.266 


9.13 




1.285 


11.1 


1.270 


9.55 




1.290 


10.5 


1.277 


11.2 




1.300 


10.3 



B-fit, s = -2.63. 



(B.2) 



The c parameters are about a factor of three smaller than for the A-fit, implying also smaller 
scaling factors A. Using the above values of s and c in a B-fit with one fit parameter tq to 
the scaling-sequence data discussed in section || leads to 



(iV 4 , k 2 ) (8000, 1.17) (16000, 1.21) (32000, 1.23) (64000, 1.26) 
r 8.87 8.96 8.67 8.65 



(B.3) 



where we have listed again the value for (JV4, K2) = (64000, 1.26). Figure 21 shows the 
A-scaled curvatures ±12/(Ar|)) for the A and B fit at N4 = 64000 as a function of «2> and 
we have also included the result of the DOB-fit in section |4.2| . Qualitatively, the switching 
from negative to positive at the phase transition is present in all fits, and in the B- and 
DOB-fits the switch is very close to k\ = 1.280(1) determined by the peak in the node 
susceptibility in p2]|. Quantitatively, the systematic differences are evidently large. 
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Figure 21: Curvatures R c = ±12/(Aro) 2 with linear interpolation. From top to bottom in k-i > 
1.28: A-fit, B-fit, DOB-fit. 



R(r) 



R(r-s) 




Figure 22: Left: curvatures R(r) from the A-fit for the K2 values of figure glj. Right: R(r — s) 
from the B-fit. A cubic polynomial was used for interpolating a(r). 



Given the parameter c = 27r 2 /v e s, the metric scale factor a(r — s) = [n(r)/c] 1 ' 5 can 



be obtained again by interpolation. The curvature R(r) following from (|4.7|) is shown in 



figure 22. As for the DOB-fits, this R(r) is not constant, although it has shallow regions 
both phases. It is evidently an observable that is very sensitive to small deviations of 
a(r) from the constant-curvature form. In particular the nice looking match of the B-fit 
to the data in figure ^0| contrasts with the non-constancy of the curves in the right plot 
of figure 22. The plots are similar to those in figure 18 in appendix [A], in which lattice 
artifacts produced by similar A- and B-fits are shown for the case of a 4D hypercubic lattice 
in flat space. Surely, there are discretization artifacts in figures ^1] and |2j, but there is 



a qualitative difference with figure 18: in the flat cubic case, the A-fit produced a fake 
positive curvature and the B-fit a fake negative curvature, whereas here the A-fit in figure 
22 shows positive as well as negative curvature, and similarly for the B-fit (right plot), with 
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an upwardly shifted R(A) compared to R(B). This supports the notion that the curvature 
plots reflect a genuine aspect of quantum spacetimes. A non-perturbative regularization in 
the continuum might may well show similar behavior. Note the small continuum distance 
r c at which the matching of continuum spacetimes to the SDT results is done: r = 8.5 in 
the matching region of the A-fit corresponds to r c = Xr ~ 4.5 and 3.8 for the crumpled and 
elongated phase; for the B-fit these are respectively r c = A(8.5 — s) ~ 3.5 and 3.8. These 
r c are just a little larger than the lattice spacing £ = \/W£ ~ 3.2. 

The results in this appendix do not change much upon relaxing some chosen constraints 
on the parameter freedom. Instead of using one value of s it can be set to the zero so found 
from the linear extrapolation of nfr) 1 ' 3 at each ki individually as in figure |I|, and also c 
can be fitted individually. This gives very similar results, with values of c differing only by 
a few percent within a phase. Of course, the individual fits improve somewhat, although 
hardly visible. However, for N± smaller than 32000 this improvement becomes important. 

C. Quartic effective potential model 

The relevant solution of the semiclassical equation can be expressed as 

-i In 

<>) = J R± + |j(-l-y±), y ± =(-Kt± + y/l + KH£) . (C.l) 

Then i? av , \ an d Xb follow from lnz± = Vw±, 

R av = r + w' + + r_w'_, r± = z±/z, (C.2) 

X = X P + Xb, (C.3) 

Xp = V [ r +(l — r +) w '+ +r_(l — r-)w'_ — 2r+r-w' + w'_] , (C4) 

Xb = r+w'l + r-w'L, (C.5) 

where the prime denotes differentiation with respect to k. Writing R± = R s ± R^, it turns 
out that Xb depends only on s± and t± (apart from n), Xp depends also on R^ but not on 
R s , and R av in addition on R s . In the Gaussian model ( 7.10| ) r_ = 1 — r + . 
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